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Intermediate-age Stars as Origin of Low Stellar velocity Dispersion 
Nuclear Rings: the case of Mrk 1157 
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ABSTRACT 

We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to map the age 
distribution of the stellar population in the inner 400pc of the Seyfert 2 galaxy Mrk 1157 
(NGC591), at a spatial resolution of 35 pc. We have performed wavelet and principal com- 
ponent analysis in the data in order to remove instrumental signatures. An old stellar popu- 
lation component (age >5Gyr) is dominant within the inner « 130pc which we attribute to 
the galaxy bulge. Beyond this region, up to the borders of the observation field, young to in- 
termediate age components (0. 1-0.7 Gyr) dominate. As for Mrk 1066, previously studied by 
us, we find a spatial correlation between this intermediate age component and a partial ring 
of low stellar velocity dispersions (en,). Low-cr, nuclear rings have been observed in other 
active galaxies and our results for Mrk 1 157 and Mrk 1066 reveal that they are formed by in- 
termediate age stars. Such age is consistent with a scenario in which the origin of the low-cr* 
rings is a past event which triggered an inflow of gas and formed stars which still keep the 
colder kinematics of the gas from which they have formed. No evidence for the presence of 
an unresolved featureless continuum and hot dust component - as found in Mrkl066 - are 
found for Mrk 1157. 

Key words: galaxies: individual (NGC591) - galaxies: individual (Mrk 1157) - galaxies: 
Seyfert - galaxies: stellar content -galaxies: near-infrared - data analysis: wavelet transform 
- data analysis: principal component analysis 



1 INTRODUCTION 

The co-existence of an active galactic 
young stars in the central region of 
wi dely known phenomenon, (e. 



200G, 



2004 



2001 ; Gonzalez-Delgado et al 



Asari et al.1 120071: iDors et al.1 120081), supporting 



nucleus (AGN) and 
Seyfert galaxies is 
Storchi-Bergmann et al. 
200 ll : Icid Fernandes et al.l 
the so- 



called AGN-St a rburst c onnection (e.g. iNo rman & Scoville 
Terlevich etafl 1 1 9901 iHeckman et al.l Il997l : iHeckmanl 



1988 



2004 



Riffel et alj|2009ch . The above studies have pointed out that the 
main difference between the stellar population (SP) of active 
and non-active galaxies is an excess of intermediate age stars in 
the former. In addition, near-infrared (NIR) SP studies have re- 
vealed that the continuum is also dominated by the co ntribution of 
intermediate - age stellar population comp onents (SPCs, Rif fel et all 
l200ll2009d . l2010d : lMartins et alj201oh . In the NIR, another com- 
ponent is commonly detected in the nucl ear spectra of Sey fert 
galaxies: unresolved hot dust emission (R iffel et al. 200 9b!ldlct 



Rodriguez- Ardila et al. 2006; Rodrigu ez-Ardila. Riffel & Pastorizal 
2005b Thus, the study of the contribution of the SPs and other corn- 
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ponents to the circumnuclear continuum of active galaxies is a fun- 
damental key in the understanding of the nature of their central 
engine. 

Very recently, the use of integral field spectroscopy with adap- 
tive optics at the Gemini North Telescope has allowed us to derive 
the contribution of distinct SPCs to the NIR spectra of active galax- 
ies. In addition, we have mapped SPCs spatial distributions and 
performed the first two-dimensional (2D) SP synthesis in the NIR 
of the nuclear region of an active galaxy (Mrk 1066, Riff el et all 
l2010bl) . A spatial correlation between the intermediate age SPC 
and a partial ring of low stellar velocity dispersions (a,) was found, 
supporting the interpretation that the low-cr, structures commonly 
observed in the inner few hundreds parsecs of Seyfert galaxies are 
due to colder reg ions with more recent star formation than the 
underlying bulge llBarbosa et alJbOOd iRiffel et alJkOQS A l2009d : 
iRiffel & Storchi-Bergmannll2010dh . Using a different method — 
modelling th e Bry equival e nt wid th, supernova rate and mass-to- 
light ratio — iDavies et al.l (2007) have quantified the star forma- 
tion history in the centre of 9 nearby Seyfert galaxies using their 
stars code. They found that the ages of the stars which contribute 
most to the NIR continuum lie in the range 10-30 Myr, pointing 
out that these ages should be considered only as "characteristic", as 
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they have not performed a proper spectral synthesis, suggesting that 
there may be simu ltaneously two or more SPs that are not coeval 
dDavies et al.ll2007l . l2006t) . 

As part of an ongoing project aimed at mapping the 2D age 
distribution of the NIR SP in the inner few hundred parsecs of 
Seyfert galaxies, we present the NIR SP synthesis of the Seyfert 2 
galaxy Mrk 1 157 (NGC 591), an early type barred spiral galaxy lo- 
cated at a distance d = 61.1 Mpc (from NASA/IPAC Extragalactic 
Database, NED), for which 1" corresponds to 296 pc at the galaxy. 
The near-IR gas excitation and kinematics as well as the stellar 
kinematics of the inner 450 pc of Mrk 1 157 are discussed in Riffel 
& Storchi-Bergmann (in prep.) using the same data used here. 

This paper is structured as follows. In Sec. |2]we describe the 
observations, data reduction and cleaning procedures. The spectral 
fitting procedures are discussed in Sec. [3] In Sec.|4]we present our 
results, which are discussed in Sec. [5] The conclusions are pre- 
sented in Sec. [6] 



2 THE DATA 

2.1 Observations and Data Reduction 

Mrk 1 157 was observed with Gemini NIFS dMcGregor et al.l2003l) 
operating with the Gemini North Adaptive Optics system ALTAIR 
in September/October 2009 under the programme GN-2009B-Q- 
27, following the standard Object-Sky-Sky-Object dither sequence, 
with off-source sky positions since the target is extended, and indi- 
vidual exposure times of 550 s. 

Two sets of observations with six on-source individual ex- 
posures were obtained at different spectral ranges: the first at the 
J-band, centred at 1.25 fim and covering the spectral region from 
1.14/jm to 1.36/im, and the second at the Ki-band, centred at 
2.3 yum, covering the spectral range from 2. 10 /mi to 2.53 pm. At 
the J-band, the J.G5603 grating and ZJ.G0601 filter were used, re- 
sulting in a spectral resolution of « 1.8 A, as obtained from the 
measurement of the full width at half maximum (FWHM) of ArXe 
arc lamp lines. The Kj-band observations were obtained using the 
KLG5607 grating and HK_G0603 filter and resulted in a spectral 
resolution of FWHM« 3.5 A. 

The data reduction was accomplished using tasks contained 
in the nifs package which is part of gemini irai[J package, as well 
as generic iraf tasks. The reduction procedure included trimming 
of the images, flat-fielding, sky subtraction, wavelength and s- 
distortion calibrations. We have also removed the telluric bands and 
flux calibrated the frames by interpolating a black body function to 
the spectrum of the telluric standard star. 

The angular resolution obtained from the FWHM of the spa- 
tial profile of the telluric standard star is 0? 1 1±0'.'02 for the J-band 
and 0'.' 12+0'.' 02 for the Kpband, corresponding to 32.6±5.9 and 
35.5+5.9 pc at the galaxy, respectively. Since the star has been ob- 
served for a shorter time than the galaxy, this value should be con- 
sidered as a lower limit for the spatial resolution. Nevertheless, the 
performance of ALTAIR for larger integration times (such those 
used for the galaxy) in previous works by our group was found to 
be similar to the one reached in the observation of the star. In the 
case of NGC4151, the spatial resolution of the galaxy is only Cf.'02 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



worse than those obtained from the F WHM of the brightness profile 
of the star ( Storchi-Bergm ann et al E009|). For each band, individ- 
ual datacubes were created at an angular sampling of O , .'05x0705 
and combined to a single datacube using the gemcombine iraf task. 
More details about observations and data reduction can be found in 
Riffel & Storchi-Bergmann (in prep.). 

2.2 Removing noise effects: The WPCA Method 

Because of the complexity of the instruments, hyperspectral data 
frequently contains instrumental signatures that cannot be removed 
with standard data reduction techniques, and which in many cases 
contaminates the science data. The Gemini NIFS is no exception. 
The wavelet principal component analysis (WPCA) technique de- 
scribed below (a full account of the technique will be addressed in 
a paper in preparation - Ferrari et al. 201 1) can separate the instru- 
mental fingerprint and, after its removal, allow the analysis of the 
"clean" data. 



2.2.1 Principal Component Analysis 

Besides the relevant information, any set of data has some amount 
of redundant information and noise. The goal of principal compo- 
nent analysis (PCA) is to find a new basis on which the data - the 
spectra from each spatial pixel - is expressed in a more meaningful 
way. The basis vectors - the principal components - are obtained 
by searching for the maximum variance in the data. There are as 
many vectors as the original variables, but usually the first few vec- 
tors contain most of the information of the data set: they provide a 
more concise description of the data. Often the physical content is 
restricted to few vectors and thus interpretation can be easier. Math- 
ematically we proceed as follows. In a basis where the variables are 
all uncorrelated (orthogonal) their covariance or correlation matrix 
is diagonal. We then find the principal components by diagonaliz- 
ing the covariance or correlation matrix of the original data. In the 
case of covariance matrix, the eigenvectors are the principal com- 
ponents and the eigenvalues are the variance associated with each 
of the eigenvectors. In the case of a hyperspectral cube, it is infor- 
mative to measure the correlation between each principal compo- 
nent and the spectra in each spatial pixel: the tomograms. For ex- 
ample, tomogram 1 is the projection (scalar product) of eigenvector 
1 and the spectra relative to each spatial pixel. For a full descrip- 
tion o n the method applied to astronomical data, see lSteiner et alj 
j2009h . 

2.2.2 Discrete Wavelet Transform 

The Discrete Wavelet Transform (DWT) consists of describing a 
signal Co in term s of a smoothed component C / and det ail wavelet 
coefficients {Wj} dMallatll999l ; IStarck & Murtaghll2006l) : 

Co = Cj + 2 Wj. (1) 

The maxi mum wavelet scale is J. In the a trous (i.e. with holes) 
algorithm dStarck & Murtagh1 l2006l) used in this work, a series of 
smoothed versions [Cj] of Co are calculated by convolving Co with 
scaled versions of a low-pass filter h. In this way, the term Cj is 
calculated by taking adjacent pixels from Cy_i that are 2 J pixels far 
apart. We proceed in this way for each of the scale j = 1 ... J. After 
this, the discrete wavelet transform is obtained from the difference 
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Figure 1. Wavelet and PCA combined: each spatial frame of the original 
datacube is first decomposed in wavelet (W) space and then the PCA analy- 
sis is performed. After identifying and correcting for unwanted effects (in- 
strumental signatures, fingerprints or noise, for example) the cube is recon- 
structed in both PCA and wavelet spaces. 

» C C • (2) 

The wavelet coefficients {Wj} now retain only details which are 
about 2 ; in size, a multiresolution transform. The maximum scale 
J is mathematically arbitrary because, by virtue of Eq. (0 and {2^, 
reconstruction is complete for any J. Physically, J is chosen so that 
the spatial frequencies of interest are lower than 2 J . 



2.2.3 The WPCA Algorithm 

The idea behind the use of the wavelet and the PCA transforms 
together is to decompose the original signal both in the wavelet 
and in the PCA space, select unwanted features and finally recon- 
struct the signal after eliminating these features. The algorithm is 
schematically shown in Figure [T] The original datacube (orange) is 
first decomposed in wavelets scales Wo, Wi, . . . , W^. Then, the PCA 
analysis is performed in each of the wavelet components, resulting 
in several eigenvectors for each scale. After removing those eigen- 
vectors which contain unwanted structures (represented by black 
dots in the figure), each of the wavelets components is restored 
with the PCA reconstruction, PCA~\ and then the resulting cor- 
rected wavelets W' Q , W[, . . . , W' k are combined to form the resulting 
corrected cube (green). 

2.2.4 Application to Mrk 1157 

Figure [2] shows an example of the WPCA transform of Mrk 1 157 
datacube. The columns are the wavelet decomposition of the orig- 
inal cube. In both FigQ] and [2] the dotted rectangle represents the 
same kind of information: in each wavelet scale a PCA analysis is 
performed. The original information is completely represented in 
the wavelet space, thus summing up all the components gives back 



the original data. The WPCA cleaning process is basically to iden- 
tify the components (in wavelet and PCA space) which contains 
mostly instrumental signatures (i.e noise), then remove them and 
reconstruct the signal. 

In the case of Mrk 1157 observations, the WPCA decompo- 
sition reveals that the Wo tomograms are dominated by the instru- 
mental signature (horizontal and vertical stripes in Fig. [2}- That 
instrumental fingerprint can also be identified in the eigenvectors 
as high frequency bumps (not shown) in the beginning and in the 
end of the spectra - the instrumental signature has a spatial and 
spectral behaviour. Thus, we remove Wo and reconstruct the cube 
with W1...W4, W c . We call the attention to the fact that we have 
analysed 5 wavelet scales and 20 principal componentes (with 20 
eigenvectors and 20 tomograms; a total of 200 maps) to identify the 
instrumental signature. For sake of brevity, in Fig. |2]we only show 
the first 4 principal components, which account for >95% of the 
data variance. 

The clean datacube was used to perform the stellar population 
analysis. The spatial resolution of the J and K; bands are very sim- 
ilar, allowing us to combine both datacubes into a single datacube 
and then perform the spectral synthesis. The combination was done 
using the scombine iraf task at a spectral sampling of 2 A. The as- 
trometry was done using the peak of the continuum emission in 
both bands. We have re-binned the resulting datacube to a spatial 
sampling of 07 1 in order to reach a signal-to-noise ratio (S/N) high 
enough to obtain a reliable fit of the stellar population. The result- 
ing S/N ranges from S/N x 20 at the borders of the NIFS field to 
up to 120 at the nucleus. We used the wspectext iraf task to convert 
the individual spectra of the final datacube to 784 ascii files, used 
as input to the synthesis code. The observation field is 2'.'8x2'.'8, 
which corresponds to 830x830 pc 2 at the galaxy. 



3 SPECTRAL SYNTHESIS 

Our main goal in this section is to map in detail the NIR spec- 
tral energy distribution (SED) components of the inner 400 pc of 
Mrk 1157. For this purpose we fit simultaneously the underlying 
continuum of the J and K bands applying the method described in 
iRiffel et al.1 J2009d) and summarized below. 

In order to p erform the stellar population synthesis we use 
the starlight code Jcl d Fernandes et al. 2004, 2005a; Mateu s"et al.l 
120061 : lAsari et alj|2007l ; Tcid Fernandes et al.ll2008l) . starlight fits 
the whole underlying spectrum, excluding emission lines and spu- 
rious data, mixing computational techniques originally developed 
for semi empirical populati on synthesis with ingredients from evo- 
lutionary synthesis models Jcid Fernandes et akl l2004, 2005a|). Es- 
sentially, the code fits an observed spectrum, A , with a combi- 
nation, in different proportions, of N* single stellar populations 
(SSPs). Extinction is parametrised by the V-band extinction A v and 
modelled by starlight as due to foreground d ust. In the fits we 
use the CCM dCardelli. Clayton & Mathislll989h extinction law. In 
order to model a spectrum M A , the code solves the following equa- 
tion: 

M A = M A0 2] Xj b hA r A (g> G(v*,o-*) (3) 

where x is the population vector, whose components xj, (j = 
1,..,N*) represent the fractional contribution of each SSP in the 
base to the total synthetic flux at A . b jA is the spectrum of the 
jth SSP of the base of elements normalized at A , the reddening 
term is represented by r A = 10~ oa(Aa ~ Aao> , M A o is the synthetic flux 
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Figure 2. WPCA decomposition of the Mrk 1157 J- band datacube. Tomograms are separated in rows, from top to bottom, tomograms 1 to 4. Original data 
(orange rectangle), is wavelet decomposed (black dotted-rectangle) in 5 scales plus continuum - w;, f = ... 4 and wc - and then reconstructed without 
the first component wo (green rectangle). Colors in individual images are false to enhance details. The first 4 tomograms shown correspond to ~ 95% of 
the variance of the data cube. Figure frames have the same color (rectangles) scheme as in FigH Color may appear inverted because tomograms, as well as 
eigenvectors, are determined from the PCA except for its sign. 



at the normalisation wavelength, denotes the convolution oper- 
ator and G(v t ,cr t ) is the Gaussian distribution used to model the 
line-of-sight stellar motions, which is centred at velocity v* with 
dispersion cr*. However, note that due to the low spectral resolu- 
tion of EPS models in the NIR, the cr* values cannot be derived in 
a reliable way from the synthesis, thus, we do not use them. For 
details on cr* see Sec. [4] The final fit is carried out minimizing the 
equation: 



XT 



(4) 



where emission lines and spurious features are excluded from the 
fit by fixing w A =0- 

In Eq. [3] the most important ingredient in stellar popula- 
tion synthesi s is the base set, b uA . As default base starlight uses 
the SSPs of iBruzual & Chariot! J2003l) . However, these SSPs do 
not include the effect of the contribution of thermally pulsating 
asymptotic giant branch (TP-AGB) stars, whose contribution is en- 
hanced in the NIR and c rucial to model the stella r populations in 
this spectral region (see iRiffel et al.l l200l l2008bl . 120094 [ 2 1 ob I 
iMartins et al.ll201oT) . Thus, we update the base using the lMarastonl 
d2005h E volutionary Population Synthesis (EPS) models as de- 
scribed in lRiffel et alj (2009d). The base comprises SSPs synthetic 
spectra covering 12 ages (t =0.01,0.03, 0.05, 0.1, 0.3, 0.5, 0.7, 1, 
2, 5, 9 and 13 Gyr) and four metallicities (Z =0.02, 0.5, 1, 2 Z Q ). 
We also include black-body fu nctions for tempera tures in the range 
700-1400 K in steps of 100 K tRiffel et alj2009dl) and a power-law 



3 ) in order to account for possible contributions from 



dust emission (BB) and from a featurel ess continuum (F C), respec- 
tively, at the nucleus (e.g. ICid Fernandes et aT, I 111). The same 
spectral base was used to map the age dis tribution of the ste llar 
population in the inner 300 pc of Mrk 1066 jRifTel et alj2010bt) . 



4 RESULTS 

In Fig. [3] (top-left) we show an optical image of Mrk 1157 
taken with the Wide Field and Planetary Camera 2 (WFPC2) 
at the Hubble Space Telescop e (HST) through the filter F606W 
dMalkan. Goriian & Tamil 19981) . The 2.17 yum NIFS datacube con- 
tinuum image is also shown in Fig. [3] (top-right). In order to illus- 
trate the accuracy of our fits, we show in Fig. [3] (bottom) sample 
spectra, obtained within O'.'lxO'.'l apertures for four distinct posi- 
tions: the nucleus (position N marked at top-left panel of Fig. 
and U.'4 west (position A), 075 south-east (position B) and 078 
north-east (position C) of the nucleus. The synthetic spectra were 
overploted on the data as dotted lines. As can be observed in this 
figure, the modelling of the stellar population reproduces very well 
the co ntinuum/absorption spectra at all positions of Mrk 1 157. Fol- 
lowing lRiffel et alj j2010bl) the observed and synthetic spectra were 
normalized at 2.1 2 fim, a region free of emission/absorption lines 
dRiffel et alj|2008bl) . 

No contribution of the FC and/or BB components was neces- 
sary in order to reproduce the the nuclear continu um in Mrk 1 1 57. 
Similar result was obtained by Riffel et al. (2009d) when fitting the 
nuclear integrated spectrum of this source within a larger aperture 
(r=:200pc) over the 0.8 — 2.4yum spectral region. We point out that 
the synthetic spectra was reddened using the A v derived from the 
synthesis (Fig.|7J. 

Following IRiffel et alj ( 1201 Obi see also Riffel et al., 2009 and 
Cid-Fernandes et al. 2004), we have binned the contribution of the 
SPCs (xj) into a reduced population vector with four age ranges: 
young (Xy\ t < 100 Myr); young-intermediate (*„: 0.3 < t < 
0.7 Gyr), intermediate-old (x io : 1 < t < 2 Gyr) and old (x a : 

5 < t < 13 Gyr). In Fig. [5] we show the spatial distribution of the 
percent flux contribution at 2.12yum of the stars in each x: while 
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Figure 3. Top-left panel: HST WFPC2 continuum image of Mrk 1157 obtained through the filter F606W (Malk an. Goriian & Tamil 19981) . Top-right panel: 
2.17 fim continuum image obtained from the NIFS datacube. Bottom panels show typical spectra obtained within an (X^SxC/.^S aperture for the nucleus and 
for a location at (X'4 north-west from it (position A). The box in the HST image shows the NIFS field of view. Bottom panels show typical spectra obtained 
within an O.' 1 xff.' 1 aperture and the resulting fit of the stellar population for the positions marked at the top-right panel. 
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Figure 4. Comparison of the overal l shape of our NIFS spectra with the 
SpeX c ross dispersed spectrum of Riffel, Rodriguez-Ardila, & P astorizal 
l2006al) . 



in regions farther than r~ 0.8(240 pc) from the nucleus the con- 
tribution of the young-intermediate SPCs reaches values of up to 
100%, closer to the nucleus the contribution of this component is 
negligible. Within this region the stellar population is dominated 
by the old component(»50%), closely surrounded (r < 01' 4) by an 
intermediate-old population (also «50%). No signs of young stellar 
populations were detected close to the nucleus; however, a signifi- 
cant contribution (25-60%) is found in regions farther than r ~0.'6. 

The light-fraction SPC contributions depend of the normal- 
ization wavelength and thus the comparison wi th results from othe r 
spectral regions should be done with caution dRiffel et alj|2010ch . 
However, a physical parameter which does not depend on the nor- 
malization point used in the fit is the stellar mass. The mass-fraction 
of each population vector components is show in Fig.[5](young: m y , 
young-intermediate: m yi , intermediate-old: m io and old population: 
m ). The maps of the mass-weighted SPC follow a similar distribu- 
tion to the light-weighted ones. However, in the former the contri- 
bution of the older ages is enhanced, particularly within 074 from 
the nucleus. 

Besides the SPC distributions, the starlight outputs the aver- 
age reddening of the stellar populations (Fig.|7). The highest values 
of E(B -V) = 0.7 (we used A v =3. 1E(B-V)) are reached at the nu- 
cleus up to within ~0.'2 from it. 

The goodness of the fit is measured in starlight by the percent 
mean deviation: adev=\Q A -M A \/ Q A , where Q A is the observed spec - 
trum and M A is the fitted model (Cid Fernandes et al .I2004 r2005b). 
The adev map for Mrk 1 157 is shown in Fig.|7]and presents values 
adev < 5 % at most locations, indicating that the model reproduces 
very well the observed underlying spectra. 

We point out that, in order to have a robust result on stel- 
lar population fitting with starlight it is important to have a 
reliable flux calibration (see starlight manual). In other words 
the fit depends on the overall shape of the observed spec- 



trum from the J to the K-bmid, and our NIFS observations 
misses the /f-band. In order to verify the reliability of our flux 
calibration we have extracted a spectrum from our datacubes 
matching the apert ure and position angle of our SpeX cross - 
dispersed spectrum (Riffel, Rodriguez-Ardila, & Pastoriza 2006a). 
Note that the SpeX data were taken in the cross-dispersed 
mode, and thus are free from seeing and aperture effects (see 
Riffel, Rodriguez-Ard ila! & Pastorizal2006al for details). The com- 
parison between both spectra is shown in Fig. |4] It is clear that our 
NIFS spectra have a reliable flux calibration. In fact the difference 
between the NIFS and SpeX spectra is lower than 3% in the J band. 
In addiditon, we have performed a simulation varying the flux of 
the / band by 20% (from 90% up to 110%) and the difference to 
the population vector components is, in mean, <5% indicating that 
our results are robust. 



5 DISCUSSION 

In general, our results are si milar to those ob tained in the 2D map- 
ping of the SP of Mrk 1066 dRiffel et alj2010bl) as well as for previ- 
ous NI R studies using single aperture nuclear spectra dRiffel et"al] 
2009d). By mapping the stellar population in 2D we can analyse 
the spatial variations of the SPCs in the inner few hundred par- 
sees of the galaxy. Further, the significance of these variations 
is enhanced by the comparison with the cr, map presented in 
Fig. [6] This map was obtained by fitting the K-band CO absorp- 
tio n band-heads with the penaliz ed Pixel Fitting (pPXF) method 
of ICap pellari & EmsellernT d2004l) using as stellar template spec- 
tra those of the Gemini library of late spectral type stars observed 
with t he Gemini Near-Infrared Spectrograph ( GNIRS) IFU and 
NIFS dWinge. Riffel & Storchi-Bergmannll2009h . More details on 
the stellar kinematics of the central region of Mrk 1157 can be 
found in Riffel & Storchi-Bergmann(i'n prep.). The cr, map shows a 
partial ring of low-cr, values (» 50-60 km s~') surrounding the nu- 
cleus at w 0/8 (230 pc) from it, immersed in higher cr, values of the 
bulge stars (» lOOkms -1 ). 

Such rings are commonly observed in the central region 
of active galaxies and are due to kinematically colder regions 
with younger stars than the underlying bulge ( Barbosa et al. 
1 2004 (Deo Crenshaw & Kraemej|2006l : ISimoes Lopes et alfeOOTt 
iRiffel et al l 12008 j. l2009al) . The comparison between Mrk 1157 
stellar population synthesis maps (light- and mass-weighted) and 
the cr» map shows that the low cr, ring is spatially correlated with 
the young-intermediate age SPC, while the highest cr,s are asso- 
ciated with th e old component. I nterestingly, very similar results 
were found bv lRiffel etal] d2010bh for Mrk 1066. Thus, the results 
found for Mrk 1157 and Mrk 1066 support the use of low stellar 
velocity dispersion as a tracer of young-intermediate age stars in 
the galaxy centre, thus confirming the interpretation of the above 
studies. 

The most compressed form to represent the SPCs mixture 
of a galaxy is by i ts mea n stellar age, which was defined by 
Cid Fernandes et al. (2005b) in two ways: the first is weighted by 
light fraction, 



(logOz. = ^x ; log; ; -, 



7=1 



and the second, weighted by the stellar mass, 



<logf*> M = ^mjlogtj. 



(5) 



(6) 
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Figure 6. Stellar velocity dispersion map from Riffel & Storchi-Bergmann 
(in prep.). 

Being both definitions limited by the age range used in our base, 
the former is more representative of the younger ag es while the 
latter is enhanced by the old SPC dRiffel et al] |2009d). In order to 
compare our results with those derived using integrated spectra we 
have calculated the mean ages for the central region of Mrk 1157. 
These mean ages in light- and mass-fractions are (logfy,) L = 8.48 
(log?,t) M = 8.65, respectively. Note that these values were obtained 
taking the mean value over all the spectra in the map (e.g. summing 
(logf*);, values over the map and dividing the result by the number 

of spectra in the map). 

Besides our study of Mrk 1066 dRiffel et alJboiObl) . to date, 
the only previous 2D stellar population studies of act ive galaxies 
in the NIR are those from the group of R. I. Davies. iDavies et al.l 



d2007h investigated the circumnuclear star formation in nine Seyfert 
galaxies using the NIR adaptive optics integral field spectrograph 
sinfoni at the Very Large Telescope (VLT). Based on measure- 
ments of the Bry emission-line equivalent width, supernova rate 
and mass-to-light ratio, these authors found circumnuclear disks 
of typical diameters of tens of pc with a "characteristic" age in 
the range 10-300 Myr. Such a "characteristic" age can be associ- 
ated with our mean age (~ 300 Myr), thus our result for Mrk 1 157 
is in good agreement with those found bv lDavies et alj d2007l) for 
their sample ( which does not inc lude Mrk 1157). Nevertheless, as 
pointed out bv lRiffel etal] d2010bl) the methodology adopted here 
allowed us not only to obtain a "characteristic" age, but to also map 
the spatial distribution of the SPCs of different ages in the central 
region of a Seyfert galaxy, on the basis of NIR integral-field spec- 
troscop y. Moreover, our res ults are in reasonable agreement with 
those o f lRiffeletalJd2009dh where the authors use integrated long- 
slit spectroscopy (considering the different apertures). 

In further support to the results of the synthesis, we found 
that the average reddening map derived for the stellar population 
(Fig.[7J is in close agreement with the one derived for the narrow- 
line region using emission-line ratios (Riffel & Storchi-Bergmann, 
in prep.), presenting a similar struc t ure and a very similar average 
value to that found by iRiffel et aT] d2009dl . E(B-V)«0.63) for the 
SP, using NIR integrated spectra. 



6 CONCLUSIONS 

As part of an ongoing project aimed at mapping the age distribution 
of the NIR SP in the inner few hundred parsecs of Seyfert galax- 
ies, we present the spectral synthesis for the nuclear region of the 
Seyfert 2 galaxy Mrk 1157 (NGC591) within the inner ~ 400 pc 
at spatial resolution of as 35 pc. We have used a cleaning method 
which allowed us to remove the noise effects and redundant data 
from the ID extractions, thus allowing a more accurate analysis of 
the SPCs. Using the noise free cube we have mapped the distribu- 
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Figure 7. Reddening map and adev map (percent mean deviation from the spectral fit). 





tion of stellar population components of different ages and of their 
average reddening. The main conclusions of this work are: 

• The cleaning technique presented here has allowed us to per- 
form the spectral synthesis in noisy spectra (e.g. borders of the dat- 
acube) in a more accurate way. 

• The age of the dominant stellar population presents spatial 
variations: the flux and mass contributions within the inner as 160 pc 
are dominated by old stars (t >5 Gyr), while intermediate age stars 
(0.3 < t < 0.7 Gyr) dominate in the circumnuclear region between 
170 pc and 350pc. 

• As for Mrk 1066 dRiffel et aI1l2010bl) . we found that there is a 
spatial correlation between the distribution of the intermediate age 
component and low stellar velocity dispersion values, which delin- 
eate a partial ring around the nucleus of Mrk 1157. Similar struc- 
tures have been found around other active nuclei, and our results 
for Mrk 1157 (and Mrk 1066) reveals that these nuclear rings are 
formed by intermediate age stars. 

• No signatures of non-thermal and hot dust components are 
found in central region of Mrk 1 157. 
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